Preconditioning Techniques for the Bidomain Equations
نویسندگان
چکیده
In this work we discuss parallel preconditioning techniques for the bidomain equations, a non-linear system of partial differential equations which is widely used for describing electrical activity in cardiac tissue. We focus on the solution of the linear system associated with the elliptic part of the bidomain model, since it dominates computation, with the preconditioned conjugate gradient method. We compare different parallel preconditioning techniques, such as block incomplete LU, additive Schwarz and multigrid. The implementation is based on the PETSc library and we report results for a 16-node HP cluster. The results suggest the multigrid preconditioner is the best option for the bidomain equations.
منابع مشابه
Scalable parallel preconditioners for an open source cardiac electrophysiology simulation package
The numerical solution of the bidomain equations is a fundamental tool in the field of computational cardiac electrophysiology. The multi-scale nature of the processes involved and the increasing complexity of the tissue and ionic models of interest make the use of High Performance Computing techniques mandatory. The solution of the linear system arising from the Finite Element discretisation o...
متن کاملPreconditioning the bidomain model with almost linear complexity
The bidomain model is widely used in electro-cardiology to simulate spreading of excitation in the myocardium and electrocardiograms. It consists of a system of two parabolic reaction diffusion equations coupled with an ODE system. Its discretisation displays an ill-conditioned system matrix to be inverted at each time step: simulations based on the bidomain model therefore are associated with ...
متن کاملFast Structured AMG Preconditioning for the Bidomain Model in Electrocardiology
The electrical activity of the heart may be modeled by a non-linear system of partial differential equations known as the bidomain model. Due to the rapid variations in the electrical field, accurate simulations require a fine-scale discretization of the equations and consequently the solution of large severely ill-conditioned linear systems at each time step. Solving these systems is a major b...
متن کاملA model-based block-triangular preconditioner for the Bidomain system in electrocardiology
We introduce a preconditioner for the solution of the Bidomain system governing the propagation of action potentials in the myocardial tissue. The Bidomain model is a degenerate parabolic set of nonlinear reaction-diffusion equations. The nonlinear term describes the ion flux at the cellular level. The degenerate nature of the problem results in a severe ill conditioning of its discretization. ...
متن کاملSolving large systems arising from fractional models by preconditioned methods
This study develops and analyzes preconditioned Krylov subspace methods to solve linear systems arising from discretization of the time-independent space-fractional models. First, we apply shifted Grunwald formulas to obtain a stable finite difference approximation to fractional advection-diffusion equations. Then, we employee two preconditioned iterative methods, namely, the preconditioned gen...
متن کامل